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Measurement of as in Radiative Hadronic Events at OPAL 
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Hadronic final states with a hard isolated photon are studied using data taken at centre-of-mass energies around 
Mz with the OPAL detector at LEP. The strong coupling as is extracted by fitting event shape variables for the 
reduced centre-of-mass energies ranging from 20 GeV to 80 GeV, and the energy dependence of as is studied. 
Combining all the values using different event shape variables and energies gives: 

as{Mz) = 0.1176 ±0.0012(stat.)lo.OTg5(syst.). 



1. Introduction 

The structure of QCD predicts a decrease of 
the strong coupling constant as for high energy 
or equivalently for short distance reactions. At 
e"''e~ collider experiments the energy scale is pro- 
vided by the center-of-mass energy ^/s. The 
LEP collider covers energy ranges from 91 to 
210 GeV P2E1. The change of as at these en- 
ergy ranges is relatively small compared to lower 
center-of-mass energies and therefore it is desir- 
able to access as at a lower y/s. 
The radiation of high energy photons either 
through initial state radiation (ISR) or through 
quark bremsstrahlung (FSR) allows to access 
lower center-of-mass energies. To make sure that 
the photon does not interfere with QCD pro- 
cesses the time scale of the photon radiation has 
to be smaller than the time scale of the par- 
ton shower. With radiative multi- hadronic events 
(i.e. e+e~ qqj a measurement of as is possi- 
ble at a reduced centre-of-mass energy, ^/s'. At 
LEPl, isolated high energy photons observed in 
the detector are mainly photons originating from 
FSR. Here we report on a measurement of as 
from event shape variables of the hadronic sys- 
tem in events with observed photons in the OPAL 
experiment 0]. The OPAL detector operated at 
the LEP e+e~ collider at CERN. A more detailed 
description of the analysis can be found in |5]. 



2. Event Selection 

2.1. Hadronic Event Selection 

This study is based on a sample of 3 million 
hadronic Z° decays selected as described in jHj 
from the data accumulated between 1992 and 
1995 at centre-of-mass energy of 91.2 GeV. We 
required that the central detector and the elec- 
tromagnetic calorimeter were fully operational. 

2.2. Isolated Photon Selection 
Isolation Cuts 

A signal event is defined as an e+e~ — *■ qq event 
with an initial or final state photon with energy 
greater than 10 GeV. The selected photon candi- 
date is required to be well contained within the 
detector. In addition the candidate cluster is re- 
quired to be well isolated from other clusters and 
tracks. The angle of the candidate cluster with 
respect to the axis of any jet, af^°, has to to be 
larger than 25°. 

Likelihood Photon Selection 

Isolated photon candidates are selected by us- 
ing a likelihood ratio method with four input vari- 
ables. The first two variables are the same as de- 
scribed above, two more variables are defined to 
reduce the background from clusters arising from 
the decay of neutral hadrons. The cluster shape 
fit variable C determines the difference between 
the measured energy deposit in a electromagnetic 
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calorimeter and the expected energy deposit from 
a simulated single photon. The second new vari- 
able is a measure of the distance between the 
electromagnetic calorimeter cluster and the as- 
sociated presampler cluster. A disagreement be- 
tween data and Monte Carlo is seen for C and 
ctfet- It is related to the difficulty in predicting 
the rate of isolated neutral hadrons in the Monte 
Carlo generators, as explained in Section In 
this analysis, the rate of isolated neutral hadrons 
used in the background subtraction is estimated 
from data. 

Hadronic events with hard isolated photon 
candidates are divided into seven subsamples 
according to the photon energy for further 
analysis. Tabled shows the mean values of \/s'(= 
2i?bGam\/l — E^/Ei^eam), and the number of data 
events. 



Table 1 

The selected number of events and the mean value 
of ^/sf for each ^/s/ subsample. 



[GeV] Events \fsi Mean [GeV] 



10-15 


1200 


78.0 ± 1.7 


15-20 


764 


71.7± 1.8 


20-25 


511 


65.0 ±2.1 


25-30 


418 


57.6 ±6.7 


30-35 


383 


48.9 ±2.6 


10-15 


303 


38.5 ±3.4 


10-15 


248 


24.3 ±5.3 



2.3. Background Estimation 

As mentioned in (7|, the JETSET Monte Carlo 
fails to reproduce the observed rate of isolated 
electromagnetic clusters, both from isolated pho- 
tons and from isolated 7r°'s. Isolated neutral 
hadrons are the dominant source of background 
for this analysis, and their rate has been esti- 
mated directly from data using the following two 
methods. 

The observed likelihood distributions in the data 
in bins of photon energy were fitted with a linear 
combination of the Monte Carlo distributions 
for signal and background events which pass the 



isolation cuts and likelihood preselection require- 
ments. The overall normalization of the Monte 
Carlo distribution is fixed to the number of data 
events. The fit is a binned maximum likelihood 
with the fraction of background events as a free 
parameter. 

Alternatively the fraction of background from 
isolated neutral hadrons was estimated from 
the rates of isolated charged hadrons. When 
isospin symmetry is assumed, the rates of neutral 
pions, neutral kaons and neutrons can be esti- 
mated from the rates of charged pions, charged 
kaons and protons, respectively. According to 
JETSET tuned with OPAL data, the deviation 
due to isospin violating decay is predicted to be 
10% for pions and 5% for kaons and protons. 
This is assigned as a systematic uncertainty for 
the isolated tracks method and combined with 
the statistical uncertainty. Consistent results are 
obtained from the two methods within the errors. 

The non-radiative multihadron background 
varies between 15.7 ± 3.5% and 4.4 ± 0.8%. The 
contamination from tt background is about 1%. 
The contribution of two photon processes is less 
than 0.01% in all subsamples. 

3. Global Event Shape Variables 

The determination of as is based on measure- 
ments of event shape variables, which are calcu- 
lated from particles with momenta in an event. 
In this analysis the following event shape vari- 
ables are used: Thrust T .8. , Heavy Jet Mass 
Mh \ and the Jet Broadening variables Bt and 
Bw &. 

The event shape variables are calculated from 
tracks and electromagnetic clusters excluding 
the isolated photon candidate. Contribution of 
electromagnetic clusters originating from charged 
particles are removed by the method described 
in [U. 

Since these variables are defined in the centre-of- 
mass frame of the colliding beams, the hadronic 
system is boosted back into the centre-of-mass 
frame of the hadrons. The Lorentz boost is deter- 
mined from the energy and angle of the photon 
candidate. The contribution from non-radiative 
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hadronic events is removed statistically by sub- 
tracting the Monte Carlo distribution scaled by 
the fraction of non-radiative hadronic events and 
r pair events listed in Table ^ The effects of 
the experimental resolution and acceptance are 
unfolded using Monte Carlo samples with full 
detector simulation (detector correction). We 
refer to the distributions after applying these 
corrections as data corrected to the hadron level. 



4. Measurement of as from Event Shape 
Distribution 

The measurement of ag is performed by fit- 
ting perturbative QCD predictions to the event 
shape distributions corrected to the hadron level 
for (1 - T), Mh US, Bt and Bw The 
0(Q;g) and NLLA calculations are combined with 
the ln(i?) matching scheme. The hadronization 
correction is applied to the cumulative theoret- 
ical calculation to conserve normalization as in 
our previous analysis at centre-of-mass energies of 
130 GeV and 136 GeV and above :2 3 . JETSET, 
HERWIG and ARIADNE are used for this hadroniza- 
tion correction and JETSET is chosen for the cen- 
tral results. 

The fit to the event shape variables uses a least 
method with as{Q) treated as a free param- 
eter. When the total number of events is small, 
the differences in statistical error between bins 
with a larger or smaller number of events than 
the theoretical prediction is not negligible. To 
deal with fluctuations of the error, the value of 
the fitted theoretical distribution is used instead 
of the number of events in each bin of the data 
distribution. The statistical uncertainty is esti- 
mated from fit results derived from 100 Monte 
Carlo subsamples with same number of events as 
selected data events. The background subtrac- 
tion, detector and hadronization corrections are 
required to be small and uniform in the region 
used for the fit. 

4.1. Systematic Uncertainties 

The total systematic uncertainties where es- 
timated by adding experimental, hadronization 
and theoretical uncertainties in quadrature. The 



largest contributions originate from the choice of 
the hadronization model and the variation of the 
renormalization scale. 

4.2. Combination of as Results 

The values of ag obtained by fits to event shape 
variables at each energy are used to study the en- 
ergy dependence of as, and to obtain an overall 
combined result for as(Afz). 
Values of as for each event shape variable and 
for all event shape variables combined are fitted 
to the solution of the renormalization group equa- 
tion at NNLO described in TT. is treated 

A'/ S 

as a free parameter in the fit. When combining 
all event shape variables, the statistical correla- 
tions between results from different variables are 
taken into account. The systematic uncertainties 

(5) 

for A-— 1- from each variable are obtained by the 

MS , 

procedure described in Section \4~1\ The value of 

(5) 

A—— derived from the values of as combining the 

MS . , , . 

event shape variables is 

AgL = 0.2027±0.0141(stat.)iO:Ji30(syst.) GeV(l) 

with xVd.o.f = 33.1/6. The system- 
atic uncertainty has contributions from exper- 
imental effects (±0.0417), hadronization effects 
(±0.0668) and variations of renormalization scale 
(+0.0810,-0.0512). This value and its total er- 
rors correspond to as(A/z) = O.lUbtomol in 
NNLO. QCD prediction of as with A^ obtained 
by the fitting is shown in Figure ^ The 
value is dominated by the subsample with Vs' 
of 38.5GeV, the othe subsamples agree well with 
the prediction. 

All values of as are propagated to the energy 
scale of Mz using the equation for the evolution 
of as quoted in |14| . The combined values and 
values for individual event shape variables for all 
y/s' subsamples are combined into one value for 
all systematic variations. Since there is no statis- 
tical correlation between values for different V7' 
subsamples, a weighted mean is calculated using 
statistical uncertainties only. The systematic un- 
certainties on the combined values are obtained 
by the procedure described in Section ITTI These 
other values are consistent with each other and 
in good agreement with the present result. As 
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Vs'[GeV] 



Figure 1. Energy dependence of ag for all sub- 
samples. 



pointed out in the analysis with non-radiative 
events at LEPl [T] and LEP2 0, the as value 
obtained by fitting B\y is lower than for the other 
three variables The combined value of as{Mz) for 
all y/V subsamples and all event shape variables 
is 

a^(Mz) = 0.1176±0.0012(stat.)±o.oo85(syst.)-(2) 

The systematic uncertainty has contributions 
from experimental effects (±0.0034), hadroniza- 
tion effects (±0.0061) and variations of renormal- 
ization scale (±0.0062,-0.0049). This value is 
consistent with the result from the analysis us- 
ing non-radiative events in LEPl data, as{Mz) = 
0.120 ± 0.006. More event shape variables are fit- 
ted for the non-radiative events. If the variables 
are restricted to the set used in the present anal- 
ysis, the combined value from the non-radiative 
events is as(Mz) = 0.1155t°;™^J. 

5. Summary 

The strong coupling as has been measured 
at reduced centre-of-mass energies, \/s', ranging 
from 20 GeV to 80 GeV using the event shape of 
the hadronic system in radiative hadronic events. 
Four event shape variables, 1 — T, Mh, Bt and 
Bw of the hadronic system boosted into the 



centre-of-mass frame are fitted with 0(ag) and 
NLLA QCD predictions and values of as are ob- 
tained for several values of v^. We obtained the 
fundamental constant of QCD, ^j^, from the en- 
ergy dependence of these values of^ as. Values at 
each 's/s' are evolved to /i Mz and combined for 
each event shape variable. The combined value 
from all event shape variables and Vs' values is 
a,(A/z) = 0.1176 ±0.0012(stat.)l°;™^3(gyg|- )^ 

This agrees with the previous OPAL analysis 
with non-radiative LEPl data and the world av- 
erage PDG value. Within errors, QCD is consis- 
tent with our data sample of events with isolated 
FSR. 
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